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Careful first-principles density functional calculations reveal the importance of hexagonal versus 
cubic stacking of closed packed planes of Pd as far as local magnetic properties are concerned. We 
find that, contrary to the stable face centered cubic phase, which is paramagnetic, the hexagonal 
close-packed phase of Pd is ferromagnetic with a magnetic moment of 0.35 /^a/atom. Our results 
show that two-dimensional defects with local hep stacking, like twin boundaries and stacking faults, 
in the otherwise fee Pd structure, increase the magnetic susceptibility. The (111) surface also 
increases the magnetic susceptibility and it becomes ferromagnetic in combination with an individual 
stacking fault or twin boundary close to it. On the contrary, we find that the (100) surface decreases 
the tendency to ferromagnetism. The results are consistent with the magnetic moment recently 
observed in small Pd nanoparticles, with a large surface area and a high concentration of two- 
dimensional stacking defects. 

PACS numbers: 73.22.-f, 75.75.+a,75.50.Cc 



I. INTRODUCTION 

Despite their narrow d-bands and high densities of 
states (DOS) at the Fermi level, which favor magnetism, 
there are only three ferromagnetic transition metals in 
nature. Palladium is paramagnetic in its stable fee struc- 
ture, but with a very high magnetic susceptibility. Sev- 
eral calculations have shown that its Fermi level lies just 
above a large peak in the DOS, at the top of the d-bands. 
The DOS at the Fermi-level, of ~ 1.1 states per spin, eV 
and atom (Chen et alX and references therein) is almost 
high enough, but not quite so, to fulfill the Stoner crite- 
rion for itinerant magnetism, since the Stoner exchange 
parameter is ~ 0.73 eV2£. The calculations have also 
shown that Pd, in its fee crystal structure, becomes fer- 
romagnetic by increasing the lattice constant by just a 
few percent. All these results have stressed the subtle 
balance of magnetism in Pd. 

Therefore, we suggest that variations in the atomic ar- 
rangement can induce changes in the density of states at 
the Fermi level that, in turn, can induce magnetism. In 
this direction, it has been proposed that monatomic Pd 
nanowires are ferromagnetic, either in their energetic 5,6 
or thermodynamic 7 equilibrium length. Also, recent ex- 
perimental results 8 indicate that fee Pd nanoparticles 
with stacking faults and twin boundaries present ferro- 
magnetism. Other experiments^ on small Pd nanoparti- 
cles also indicate the existence of a hysteresis loop which, 
in this case, is interpreted as due to a non zero magnetic 
moment at the surface atoms. 

Stacking faults and twin boundaries are very common 



two-dimensional defects in face centered cubic (fee) met- 
als as well as in diamond and wurtzite-structure semicon- 
ductors. Their abundance is partly due to their low en- 
ergy of formation, since they preserve the local geometry 
and close packing. Their importance in the mechanical 
properties, like hardness and brittleness are recognized 
since many years. Recently, the role played by the stack- 
ing of (111) layers has been emphasized in connection 
with magnetic orientation and magnetic ordering of thin 
layers and superlattices of mag netic metalo 10 : 11 : 12 . The 
growth of Co on top of the Cu(lll) surface has been 
extensively studied experimentally, and the correlation 
between stacking pattern and magnetic properties has 
been reasonably well established. 

At this stage it seems worth studying the effect of 
the stacking sequences and of the surfaces on the elec- 
tronic and magnetic properties of Pd. With these ideas 
in mind, we have calculated total energies and magnetic 
moments of Pd in both the fee and hexagonal close- 
packed (hep) phases, as well as in surfaces and near two- 
dimensional stacking defects like intrinsic and extrinsic 
stacking faults, and twin boundaries. 



II. METHODOLOGY 

All our calculations are performed within density 
functional theorjii^ (DFT), using either the local den- 
sity approximation^ (LDA) or the generalized gradi- 
ent approximation 15 (GGA) to exchange and correla- 
tion. Most of the calculations were obtained with the 



2 



SIESTAiSiii method, which uses a basis of numerical 
atomic orbitals— and separable 1 ^ norm conserving pseu- 
dopotentials 20 with partial core corrections 21 . To gen- 
erate the pseudopotentials and basis orbitals, we use a 
Pd configuration 4d 9 5s 1 , since we have checked that it 
leads to better transferability and bulk properties than 
the ground state configuration 4g? 10 5s°. After several 
tests we have found satisfactory the standard double-£ 
basis with polarization orbitals (DZP) which has been 
used throughout this work. The convergence of other 
precision parameters was carefully checked. The range of 
the atomic basis orbitals was obtained using an energy 
shiftii of 50 meV. The real space integration grid had a 
cut-off of 500 Ryd, while around 9000 k points/atom -1 
were used in the Brillouin zone sampling. To accelerate 
the selfconsistency convergence, a broadening of the en- 
ergy levels was performed using the method of Methfessel 
and Paxton^ 2 " which is very suitable for systems with a 
large variation of the density of states at the vicinity of 
the Fermi level. 

It is necessary to mention that most of the energy dif- 
ferences between paramagnetic and ferromagnetic solu- 
tions in Pd structures are extremely small, what requires 
a very high convergence in all precision parameters and 
tolerances, and specially in the number of k points. It 
must be recognized, however, that the basic DFT uncer- 
tainty is probably larger than those energy differences, so 
that it is not really possible to determine reliably whether 
a particular defect or structure is para or ferromagnetic. 
Still, we think that it is possible to find reliably the rela- 
tive tendency towards magnetism of different structures. 
In particular, E(M) curves, of total energy versus total 
magnetic moment, provide an excellent tool to study the 
tendency to magnetism in different systems: indepen- 
dently of whether the systems are para or ferromagnetic, 
one may determine if a defect or a surface has a smaller 
or a larger tendency to magnetism than the bulk, de- 
pending on which of the two E(M) curves is higher. In 
addition, the selfconsistent convergengy is considerably 
faster at constant magnetic moment, so that it is possi- 
ble to determine reliably (for a given functional) relative 
energies as small as a fraction of an meV. 



III. BULK CRYSTALS 

The first mandatory system to be considered is the per- 
fect bulk crystal in its experimentally-stable fee phase. In 
principle, the GGA functional goes beyond the LDA and 
it is generally considered to be more accurate and reli- 
able. However, in the case of Pd, we find that the GGA 
gives a lattice constant of 3.99 A, 2.5% larger than experi- 
ment, and a ferromagnetic ground state with a magnetic 
moment of 0.4/is/atom and an energy 4.5 meV/atom 
below the paramagnetic phase (Figure 1). Since this 
energy difference is very small, and to discard that the 
ferromagnetic phase is favoured by the pseudopotential 
or basis set used, we have reproduced^ this result using 
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FIG. 1: Magnetic moment versus lattice constant for fee Pd 
obtained with the SIESTA code and the LDA and GGA func- 
tional. The insets indicate the variation of the total energy 
with magnetic moment at the corresponding equilibrium lat- 
tice constants, indicated by arrows. 



two other DFT methods: the pseudopotential plane wave 
code VASP— and the all-electron augmented plane wave 
method WIEN24, both with the same GGA functional^. 
Previous GGA calculations by Singh and Ashkenazi2£, 
using a different functional, also found a lattice constant 
larger than the experimental one, but did not address the 
subtle ferromagnetic-paramagnetic transition with suffi- 
cient detail. 

On the other hand, within the LDA we find a paramag- 
netic ground state and a lattice constant of 3.89 A, both 
in agreement with the experimental values and with pre- 
vious LDA calculations of Moruzzi and Marcus 4 . These 
results suggest that the GGA is not necessarily more re- 
liable to study magnetism in Pd. In fact, it is not feasible 
to study the possible existence of magnetic defects in fee 
Pd when the bulk result is already ferromagnetic. There- 
fore, we have chosen the spin-dependent LDA to perform 
most of the remaining calculations in this work, although 
GGA results have been obtained also as a check in some 
cases. 

The magnetic susceptibility \ °f the magnetic moment 
M to an external magnetic field H is 

_ fdM\ _ ( d 2 Ey 1 

Thus, the flatness of the E(M) curve in Figure 2] implies 
a very large susceptibility, even within the LDA. To com- 
pare with the other transition metals in the same column 
of the Periodic Table we have calculated the variation of 
energy with magnetic moment for Ni and Pt. The results 
are shown in Figure |21 We immediately observe a clear 
ferromagnetic and paramagnetic behavior of Ni and Pt 
respectively, while Pd, as indicated above, is paramag- 
netic but very close to the ferromagnetic transition. 

In order to study how magnetism depends on the local 
geometry and stacking of atoms, we have first consid- 
ered the hep structure as compared to the fee one. It is 
known that the breaking of the cubic symmetry in stack- 
ing faults, while keeping the number of nearest-neighbor 
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FIG. 2: Total energy versus magnetic moment for three tran- 
sition metals in the same column of the Periodoc Table, cal- 
culated in the LDA. Note the flatness of the Pd curve, re- 
sponsible of the large magnetic susceptibility of this metal. 



atoms and their bond lengths and angles, induces a rear- 
rangement of the d energy bands that can even give rise 
to localized electronic states 2 ^. The DOS of hep and fee 
structures, in the paramagnetic phase, shown in Figure 
EJ were calculated for the same nearest neighbor distance 
of 2.76 A and the ideal ratio c/a = (8/3) 1/2 for the hep 
structure, since we obtain that both structures are sta- 
ble at this distance, with the hep c/a ratio only 0.75% 
larger than the ideal value. Several points are worth men- 
tioning: (i) The d bandwidth is almost identical in both 
structures because they have identical nearest-neighbor 
configurations, (ii) The shapes of the DOS are very dif- 
ferent, as a consequence of the breaking of the cubic 
symmetry: in fee the second-nearest-neighbor atoms are 
staggered whereas in hep they are eclipsed and the lack 
of cubic symmetry inhibits the ti g — e g splitting of the 
d bands. (Hi) The DOS at the Fermi level is larger in 
hep than in fee (1.43 versus 1.15 states per atom per eV 
and per spin). This implies that Pd in the hep struc- 
ture satisfies the Stoner condition for ferromagnetism: 
with a Stoner exchange parameter 7=0.73 eV2&, we get 
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FIG. 3: Calculated electronic densities of states of Pd in 
the paramagnetic phase (i.e. forcing equal spin up and spin 
down populations) in the fee and hep crystal structures. The 
vertical line indicates the Fermi level. 



I x D(E F ) = 0.84 for fee and I x D(E F ) = 1.05 for hep. 

We have then calculated the total energy and the mag- 
netic order in the hep phase in the LDA. The results 
are shown in Figure We immediately observe that, 
contrary to what happens in the fee structure, there is 
a non zero magnetic moment of 0.35 /Lt^/atom at the 
equilibrium interatomic separation for the hep structure, 
<i=2.76 A. Like in the fee phase, the system is close to 
a magnetic-nonmagnetic transition but, in this case, in 
the ferromagnetic side, with the ferromagnetic hep phase 
approximately 1.7 meV lower in energy than the param- 
agnetic one. The transition from para to ferromagnetism, 
as a function of the lattice constant, is abrupt, like in the 
fee case, although we cannot assess with enough confi- 
dence whether it is first or second order^Si. 

Finally, the hep structure is 4.0 meV higher in energy 
than the fee. This ordering is in agreement with experi- 
ment, and we reproduce it also with the GGA, but it is in 
contradiction with the calculations of Huger and Osuch 28 
that reported an hep structure lower in energy, but also 
ferromagnetic. 



IV. BULK DEFECTS 

Hexagonal (111) closed packed planes of atoms can be 
stacked in different ways, giving rise to different struc- 
tures. If we label the three possible positions of the atoms 
as A, B and C, the fee stacking is ...ABCABC... and that 
of hep is ... ABABAB ... 29 -, where the layers with hexag- 
onal symmetry are undelincd. Different defects can be 
generated in the fee structure: The intrinsic and extrinsic 
stacking faults have two hexagonal layers, with stacking 
sequences ...ABCACABC... and ...ABCACBCABC... 
respectively. The twin boundary has a single hexagonal 
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FIG. 4: Calculated total energy for (a) fee and (b) hep Pd as a 
function of the nearest neighbor distance. The arrows indicate 
the minimum (equilibrium) nearest-neighbor distance. The 
insets at (a) and (b) show the total energy versus magnetic 
moment for the fee and hep structures respectively. In these 
cases the energy origin is at the corresponding minima. 



layer with ...ABCACBA... stacking. 

Previous model calculations^ at stacking faults of 
transition metals indicate an important perturbation of 
the local densities of states from that of the perfect fee 
lattice. The presence of a hexagonal stacking of layers in- 
duces localized electronic statesS^ . and an enhancement 
of the DOS at the top of the valence band. These calcu- 
lations suggest possible variations of the magnetic prop- 
erties around the extended defect. We have performed 
calculations of several packing sequences, in supercells 
containing various numbers of fee layers, to study to what 
extent different local configurations can give rise to mag- 
netic moments. The results of the calculations, both in 
the perfect crystals and in the defects, are independent 
(within less than 1%) of the initial input moment. 

Most of the calculated supercells have a finite mag- 
netic moment and some are shown in Figure|3J However, 
its variation with the supercell thickness is rather com- 
plex and nonmonotonic. We suspect that this complex 



FIG. 5: Local magnetic moment in supercells of different 
thickness, containing a single intrinsic stacking fault between 
fee layers, calculated with the LDA. 



behaviour is due to an oscillatory component of the mag- 
netic coupling between the neighbouring stacking faults, 
like that observed in superstructures of magnetic slabs 
sandwiched between nonmagnetic metals^I. Although we 
have not been able to stabilize any antiferromagnetic so- 
lution in double-size supercells, this may be due to the 
size limitations and to the difficulties of convergence. As 
expected, the magnetic moments are smaller in fee than 
in hep layers, but they are nevertheless far from negligi- 
ble. The interplay between the magnetism of the hexago- 
nal layers and the paramagnetism of the cubic ones, both 
close to a paramagnetic to ferromagnetic transition, is 
very subtle. What is important is that the hexagonal lay- 
ers have a tendency to become magnetic and they induce 
magnetic moments at the atoms in the cubic layers. This 
is due first, to the large magnetic susceptibility in fee Pd, 
where the hexagonal layers act as magnetic impurities. 
And second, because of the two-dimensional character 
of the defects, any magnetic perturbation decays very 
slowly with distance, like in a pseudo one-dimensional 
metal. In other words, there is a long-range RKKY-like 
interaction between the hep layers through the interme- 
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diate fee ones. 

Given the size limitations of our calculated supercells, 
and the nonmonotonic magnetic moment with increasing 
supercell thickness, it is not possible to conclude con- 
fidently whether isolated planar defects in fee Pd have 
a finite magnetic moment, within the LDA. However, 
Fig. El clearly shows that all the hexagonal defects have 
a larger tendency to magnetism than the bulk fee lat- 
tice, as indicated by their respective E{M) curves. This 
is hardly surprising, given the ferromagnetic character 
of hep Pd. On the other hand, all the supercells are of 
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FIG. 6: Total energy as a function of the magnetic moment 
for differents two dimensional defects of Pd, with local hexag- 
onal symmetry, compared to the bulk fee lattice. The energy 
at M = has been shifted to a common value. The defects 
are separated by 6 fee layers like in Figure 5(b). The stacking 
fault curve has a minimum of -0.2 meV at M ~ 0.09/1 b /atom, 
in agreement with Figure^, which is not noticeable at the 
figure scale. 

course strongly magnetic within the GGA, with a larger 
tendency to magnetism (lower E(M) curve) than the 
bulk fee crystal. Therefore, it is perfectly possible that 
an isolated stacking fault in fee Pd is indeed magnetic in 
nature. 



V. SURFACES 

Recent experimental results^ on Pd nanoparticles have 
been interpreted as magnetism at (100) surfaces. In prin- 
ciple, surface magnetism is plausible in general because 
the lower coordination of surface atoms favors narrower 
bands and larger densities of states. In practice, how- 
ever, surface relaxation and reconstruction may contract 
the surface bond lengths and more than compensate for 
the lower coordination. We have then calculated the 
electronic structure of finite Pd slabs in the (100) and 



(111) orientations, after carefully relaxing their geomme- 
try. The results of the calculated total energy versus 
magnetic moment are shown in Figure \7\ for the largest 
calculated thickness. This general tendency is consis- 
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FIG. 7: Total energy within the LDA versus magnetic mo- 
ment of bulk fee Pd and of 9-layer slabs with surfaces oriented 
in the (100) and (111) directions. The energy at M = has 
been shifted to a common value. 

tently obtained for sufficiently thick slabs, but the thik- 
ness dependence of the total magnetic moment, shown in 
Fig- IHi s n °t monotonic, as in the bulk supercells. 
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FIG. 8: Largest local magnetic moment, as a function of the 
number of layers, in slabs with surfaces oriented in the (100) 
and (111) directions. 
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We observe in Fig. [7| that the minimun energy is at 
zero magnetic moment and therefore the two surfaces are 
paramagnetic. Moreover, the E(M) curve of the (100) 
surface is higher than that of the bulk, showing that this 
surface is less prone to magnetism than the bulk. On 
the other hand, the (111) surface, although also param- 
agnetic within the LDA, has a larger susceptibility than 
the bulk, what makes plausible that it may be magnetic 
in nature. 

We have found that bulk planar defects, as well as 
(111) surfaces, have a larger tendency to magnetism than 
the bulk Pd crystal, which is itself on the verge of ferro- 
magnetism. Since both surfaces and defects are in high 
concentrations in nanoparticles, and both have in fact 
been proposed independently A£ as responsible of the 
observed magnetic moment of these particles, it makes 
sense to consider their combined effect. To this end, we 
have calculated the geometry, energy, and magnetism, of 
planar defects close to a (111) surface. We use a slab in 
which the opposite surface is "magnetically passivated" 
by imposing a short distance between the first two atomic 
planes, what inmediately kills their local magnetic mo- 
nent. In this way, we ensure that the possible magnetic 
moment of the slab is due to the combination of the stack- 
ing fault and a single surface (in fact, this structure pe- 
nalizes and sets a lower limit for the appearence of mag- 
netism). Still, as shown in Figure [5] we find a clear mag- 
netic moment in all the cases considered. This tendency 
is further demonstrated by their E(M) curves, presented 
in Figure 1101 for two cases, which show unambiguously 
their magnetic character. Furthermore, the total energy 
of all the slabs is lower than that of the defects in the 
bulk, plus that of the unfaulted slab, implying that the 
surface atracts the defects and that they should therefore 
be expected to appear together in nanoparticles, due to 
the high concentration of stacking faults and to the small 
space between surfaces. Notice that the extrinsic stack- 
ing fault is more stable than the intrinsic one both at the 
surface and at the bulk. 



VI. CONCLUSIONS 



DO 

a 

u 

£ 
o 

u 

a 

Ed 



0.3 
0.2 
0.1 
0.0 

0.3 
0.2 
0.1 
0.0 

0.3 
0.2 
0.1 
0.0 

0.3 
0.2 
0.1 
0.0 

0.3 
0.2 
0.1 
0.0 



[a) 





■ 






AE=-2.5 meV 


■ 

i 


i 




■ 

1 


■ 


f . . - 


A 

A 


(j 


A 

A 


B 




A D C* A 
A D O A 


■b) 












AE=-9.9 meV 


1 

i 


1 


i 


i 


■ 


W . j. j. 


- A 


B 


A 


B 


c 


A B C A - 




■ 


o 




i 




AE=-4.5 meV 


■ 

i 


i 


1 


■ 

1 


f m j. + 


- A 


B 


C 


B 


c 


A B C A - 


]d) 


^ 









AE=-11.1 meV 


■ 















i 


i 




i 




■ 

■ j. j. 


- A 


B 


C 


A 


c 


A B C A - 




O 


■ 







AE=-14. 


9meV 


■ 

i 


i 




i 


■ 

1 


"dead 

■ 

* 


ayep 
i 


C 


A 


c 


B 


c 


A B C A - 



FIG. 9: Local magnetic moments for diferents stackings of 
hexagonal and cubic planes in a slab of 9 layers, (a) Twin 
boundary, (b), (c), and (d) Intrinsic stacking faults, with the 
hep layers in differents positions, (e) Extrinsic staking fault. 
Squares and hexagons represent cubic and hexagonal layers, 
respectively. The distance between the two rightmost layers 
was fixed to a small value, in order to kill their tendency 
to magnetism and thus to simulate the bulk. The energies 
reported are the difference between the total energy of the 
slab, minus those of the defect in the bulk and of the unfaulted 
slab. The formation energies in the bulk are 33.9, 74.7 and 
73.0 meV for the twin boundary, intrinsic stacking fault, and 
extrinsic stacking fault, respectively. 



Our main conclusions can be summarized as follows: 
(i) The GGA functional gives an incorrect ferromagnetic 
ground state for the fee Pd crystal. Accordingly, all the 
defects studied are also magnetic within the GGA but, 
obviously, this does not imply that they are magnetic 
in nature. On the contrary, the simpler LDA gives the 
correct lattice constant and paramagnetic state, (ii) The 
hep phase is ferromagnetic, within both the LDA and 
GGA. In the LDA, it has an energy 1.7 meV lower than 
the hep paramagnetic state and 4.0 meV above the fee 
phase. (Hi) We cannot determine whether an isolated 
stacking fault is magnetic in the LDA, but it certainly has 
a larger magnetic susceptibility than the perfect crystal, 
and might be magnetic in nature, given the uncertainty 
between the different functionals. (iv) The free (100) 



surface is paramagnetic, with a lower susceptibility than 
the bulk crystal, (v) The (111) surface is paramagnetic 
in the LDA, but it has a larger susceptibility than the 
bulk crystal, and it might be magnetic in nature, (vi) 
Hexagonal planar defects are atracted towards a (111) 
surface, and they become clearly magnetic when close 
enough. 

Therefore, (100) surfaces are not good candidates as 
the origin of magnetism in Pd nanoparticles, as they 
had been proposed^. In contrast, planar stacking de- 
fects, (111) surfaces, and specially a combination of both, 
are plausible candidates to present permanent magnetic 
moments and to be responsible for the magnetism ob- 
served in Pd nanoparticles. Thus, our results are con- 
sistent with experiments in small Pd clusters of average 
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FIG. 10: Calculated total energy, as a function of the mag- 
netic moment, for two different two-dimentional defects close 
to a (111) surface (shown in Figure|5|(c) and (e)). 



diameter 2.4 nm, which are reported to display sponta- 
neous magnetization 8 . High-resolution transmission elec- 
tron microscopy has shown that a high percentage of the 



particles exhibit single and multiple twinning boundaries. 
In addition, the smallness of the spontaneous magnetiza- 
tion seems to indicate that only a small fraction of atoms 
hold a permanent magnetic moment and contribute to 
ferromagnetism. Other experimental results 9 on small 
Pd particles have also shown their ferromagnetic charac- 
ter. Besides, ferromagnetism can also take place in other 
non ideal structures like nanowires^. 

Magnetic anomalies observed experimentally in differ- 
ent Ni2£ and Coi£ stacking can be interpreted along the 
lines described in this work. The fact that stacking- 
faults in Pd display non negligible magnetic moments, 
not present in bulk fee crystal, opens new lines of re- 
search. The appearance of magnetism in nominal fee 
samples should be revisited in view of our results, since 
so far the possibility of magnetism around stacking faults 
has been overlooked. Also, layer growth of Pd on top of 
non magnetic substrates, including Pd itself, may pro- 
duce an interesting magnetic phenomenology in connec- 
tion with the stacking structure which, in turn, depends 
on the method used for growth. The study of the depen- 
dence of the magnetic properties of Pd grown Ag(lll) 
and Pd/Ag multilayers on the stacking sequences is un- 
der way and will be reported elsewhere. 
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